arXiv:1504.04724vl [gr-qc] 18 Apr 2015 


Newtonian-noise cancellation in full-tensor gravitational-wave 

detectors 


Jan Harms 

Universita degli studi di Urbino “Carlo Bo” and 
INFN, Sezione di Firenze, Sesto Fiorentino, F50019, Italy 

Ho Jung Paik 

Department of Physies and Astronomy, 

University of Maryland, College Park, Maryland 20742 

Abstract 

Terrestrial gravity noise, also known as Newtonian noise, produced by ambient seismic and in¬ 
frasound fields will pose one of the main sensitivity limitations in low-frequency, ground-based, 
gravitational-wave (GW) detectors. It was estimated that this noise foreground needs to be sup¬ 
pressed by about 3 ~ 5 orders of magnitude in the frequency band 10 mHz to IHz, which will 
be extremely challenging. In this article, we present a new approach that greatly facilitates can¬ 
cellation of gravity noise in full-tensor GW detectors. The method uses optimal combinations of 
tensor channels and environmental sensors such as seismometers and microphones to reduce gravity 
noise. It makes explicit use of the direction of propagation of a GW, and can therefore either be 
implemented in directional searches for GWs or in observations of known sources. We show that 
suppression of the Newtonian-noise foreground is greatly facilitated using the extra strain channels 
in full-tensor GW detectors. Only a modest number of auxiliary, high-sensitivity environmental 
sensors are required to achieve noise suppression by a few orders of magnitude. 

PACS numbers: 04.80.Nn, 95.55.Ym, 07.60.Ly, 42.62.Eh, 04.80.-y 
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I. INTRODUCTION 


The advanced generation of large-scale, laser-interferometric GW detectors LIGO [1] 
and Virgo [2] are expected to make the hrst direct detections of GWs within the next 
few years, which will open a new observational window to the Universe. The Japanese 
GW detector KAGRA is currently under construction and will join the detector network 
near the beginning of the next decade [3]. These kilometer-scale detectors are designed 
to observe GWs in the frequency band between 10 Hz and a few 1000 Hz. Upgrades of 
these detectors can potentially extend the band to lower frequencies by a few Hz [4, 5], 
but today it seems infeasible to continue developing the existing facilities into detectors 
sensitive well below 10 Hz. Projecting the state-of-the-art GW detector technology into the 
near future implemented in a detector with 10 km arm lengths, and assuming a new detector 
site favorable in terms of ambient seismic noise (and associated gravity noise), it seems 
feasible to extend the detection band down to frequencies around 3 Hz, as was the result of a 
design study for the European third-generation detector Einstein Telescope [6] . Gompletely 
new detector designs need to be considered to realize ground-based GW detectors at even 
lower frequencies [7]. These include the atom-interferometric [8], the torsion-bar [9], and the 
superconducting [10] GW detector concepts targeting signals between 10 mHz and 1 Hz. 

The low-frequency sensitivity goals set for any of the potential future ground-based detec¬ 
tors is strongly influenced by estimates of Newtonian noise (NN). If an ideal site is selected, 
which means that seismic and infrasound noise are near their global minima [11], then GW 
strain sensitivities of a few times 10“^^ Hz“^/^ can be reached down to a few Hz without 
further gravity-noise mitigation techniques. At less favorable sites, such as the existing de¬ 
tector sites, or considering lower-frequency detectors, gravity-noise mitigation is required. 
Proposed strategies can be divided into two categories; passive and active noise mitigation. 
Passive mitigation aims at suppressing sources of gravity perturbation close to a detector’s 
test masses. The detector buildings hosting the test masses act as a shield against envi¬ 
ronmental infrasound suppressing associated NN [12]. The construction of moats has been 
proposed reflecting incoming seismic surface waves as a means to reduce NN at the LIGO 
sites [13]. Recess structures around the test masses can also reduce seismically induced 
NN [14]. However, as was explained in [14], these techniques are effective only at higher 
frequencies around 10 Hz, where the mitigating structures can have dimensions similar to 
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the lengths of infrasound or seismic surface waves. Site selection is also considered a passive 
mitigation strategy. Building a GW detector underground, such as the KAGRA detector, 
greatly suppresses NN from seismic surface waves above a few Hz [15]. 

Whenever the passive strategies are not an option or resulting noise suppression is in¬ 
sufficient, active noise mitigation needs to be considered. Gommon to all active mitigation 
strategies is the usage of an array of environmental sensors with the purpose to obtain in¬ 
formation about density perturbations in the vicinity of the test masses. Implementations 
of these methods then differ in how one makes use of these data. One could actively cancel 
the density perturbations near test masses using optimal feedback control. For example, 
microphones can be controlled to produce sound that cancels the ambient sound held in¬ 
side a chosen volume [16]. A similar scheme may be possible for seismic helds. However, 
this approach cannot be effective at very low frequencies where active cancellation of den¬ 
sity huctuations must be exerted over large volumes around test masses. Another idea is 
feed-forward subtraction where an estimate of the NN obtained from environmental data is 
used to cancel the gravity-induced motion of a test mass, or similarly, the estimate can be 
subtracted from the detector’s data in a post-processing step. This method was hrst inves¬ 
tigated in detail for the case of stationary NN using Wiener hlters [17], and later also tested 
in numerical simulations of non-stationary seismic helds [18]. While in the last publication 
suppression of NN from seismic surface waves was achieved using a relatively small array 
of about 10 seismometers, it is to be expected that especially subtraction of infrasound NN 
at frequencies where it is relevant (below a few Hz) requires a large number of auxiliary 
sensors [7] . Suppression of infrasound NN below 1 Hz by orders of magnitude is considered 
an extreme challenge and potential show-stopper for low-frequency GW detectors. 

In this paper, we outline a crucial advantage of full-tensor GW detectors over conventional 
detectors that only output one component or combination of components of the gravity 
strain tensor. Full-tensor detectors measure all 5 independent components of the gravity- 
strain tensor as explained in Section H for the example of superconducting GW detectors. 
The basic idea behind the new cancellation scheme is that a suitable linear combination 
of some tensor channels should make it possible to cancel NN in the remaining tensor 
channels. This problem is studied analytically in Section HI. It will be shown that tensor 
NN cancellation still requires auxiliary environmental sensors, but the problem is greatly 
facilitated by including tensor channels. The optimal combination of tensor channels depends 
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on the direction of propagation of a GW, and therefore the method can be applied in 
directional searches of GWs or observations of known GW sources. The role of sensor 
noise is emphasized in Section IV, where Wiener hlters instead of the analytical expressions 
are introduced to hnd optimal channel combinations. In Section V, we propose a practical 
implementation of the method based on the simulated noise suppression using Wiener hlters. 


II. FULL-TENSOR GW DETECTORS 

According to general relativity, a gravitational held is characterized by a curvature tensor. 
Terrestrial laser-interferometer GW detectors measure only one oh-diagonal component by 
combining two orthogonal light cavities. A full-tensor detector could be constructed by 
measuring hve degenerate quadrupole modes of a solid sphere [19, 20]. The bandwidth 
of the detector could be widened by using a ’’split sphere,” in which six test masses are 
suspended from a central mass [21], or a ’’dual sphere,” in which a spherical shell encloses 
an inner sphere [22]. A tensor detector is equally sensitive to GWs coming from any direction 
with any polarization and is thus capable of resolving the source direction and polarization. 

One could construct a low-frequency (0.01 Hz to 10 Hz) tensor GW detector by using 
six ’’almost free” test masses [10]. Figure 1 shows the test mass conhguration of such a 
detector. Six superconducting test masses, each with three linear degrees of freedom, are 
levitated over three orthogonal mounting tubes. The test masses are made of niobium (Nb) 
in the shape of a rectangular shell. Superconducting levitation/alignment coils and sensing 
capacitors (not shown) are located in the gap between the test masses and the mounting 
tubes, as well as on the outer surfaces of the test masses. The along-axis motions of the two 
test masses on each coordinate axis are differenced to measure a diagonal component of the 
wave: 

hii{u:) = ^{x+ii{u:) - X-ii{u:)), ( 1 ) 

where x±ij{uj) is the displacement amplitude of the test mass on the ±z axis along the 
j-th axis and L is the separation between the test masses on each axis. The cross-axis 
(rotational) motions of the four test masses on each coordinate plane are differenced to 
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FIG. 1. Test mass configuration for the low-frequency superconducting tensor GW detector. Mo¬ 
tions of six magnetically levitated test masses are combined to measure all six components of the 
curvature tensor. 


measure an off-diagonal component of the wave: 

- {x_ji{u) - x+ji(cn))], i ^ j. 


( 2 ) 


In addition to measuring the six strain signals, the detector will measure the three linear 
and three angular platform acceleration (plus gravity) signals by summing the along-axis 
and cross-axis test mass motions: 

aii{u:) = — —{x+ii{u:) + X-ii{u:)) 

up' 

= - Y [(^+p(^) - 

+ - x+jpuj))] ppj. 
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These common-mode (CM) acceleration signals are used to remove the residual sensitivity 
of the GW detector to the platform accelerations [10]. 

Since test mass motion is measured with respect to the sensing circuit elements mounted 
on the platform, this detector requires a rigid platform with mode frequencies above the 
signal bandwidth. To reduce its thermal noise, the platform itself may need to be cooled 
to 77 K or lower. To alleviate excessive demand on cryogenics, the platform must not 
be too heavy while it is rigid enough, with all the resonance frequencies above 10 Hz. The 
design details of this low-frequency tensor detector, called SOGRO (Superconducting Omni¬ 
directional Gravitational Radiation Observatory), will be published elsewhere (Paik et ah, 
in preparation). 


III. NEWTONIAN NOISE FROM INFRASOUND AND SEISMIC SURFACE 
FIELDS 

In this section, we present the analytical relations between NN contributions to different 
channels of the full-tensor GW detector. We consider the two cases of NN from seismic 
surface waves and infrasound, which are considered the dominant contributions to terrestrial 
gravity noise below 1 Hz. Rayleigh waves are the only surface waves producing gravity 
perturbations. The perturbation of the gravity potential above the surface produced by 
Rayleigh waves is given by [7] 

= -27r^^^^^(a;)e"^V'''^. (4) 

Here, q = {x,y), k = A;(cos(a),sin(Q!)), G denotes the gravitational constant, po the mean 
mass density of the ground, ^(cu) the amplitude of vertical surface displacement, and 7 ~ 
0.8 a numerical factor characteristic for fundamental Rayleigh waves that depends on the 
ground’s Poisson’s ratio. We will assume here that the horizontal wavenumber k obeys the 
linear dispersion relation k = cu/cRf, but this is not important for the method and only 
simplihes the equations. While the perturbed gravity potential underground has additional 
contributions, for example, from the displacement of cavity walls, hnal results presented 
in this section are independent of this as long as the depth of the GW detector is not 
a signihcant fraction of the Rayleigh-wave length. So the choice of considering a surface 
detector for the Rayleigh NN calculation is just to simplify some equations. 
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The gravity gradient tensor is defined as 


g(f,a;) = -V O V5(/)Rf(r,a;). 


(5) 


The gravity-gradient tensor can be identified with the second time derivative of gravity 
strain, g = h. This equivalence holds at low frequencies and for ground-based detectors 
where the effect of a GW can effectively be described as a tidal force acting on test masses. 
The response of low-frequency GW detectors to NN is described by this tensor since the 
distance between test masses is much smaller than the length scale of variations in the 
gravity field. In other words, the expression for strain NN in large-scale GW detectors, 
—V50Rf/L, with L being the distance between test masses, is approximately given by the 
second spatial derivative in Equation (5). Gonsequently, NN in low-frequency detectors is 
independent of L. 

For Rayleigh waves, we have 


gRf(h'= 0,a;;a) = 27i'yGpo^{u)k 

^ cos^(a) cos(a)sin(a) 
cos(q!) sin(a) sin^(a;) 

, —icos(a) —isin(a) 


—icos(a)^ 
—isin(a) 


-1 


/ 


( 6 ) 


An arbitrary Rayleigh-wave field can be written as a sum over many individual waves. 
Parameterizing the direction of propagation of a GW by angular spherical coordinates 0, 0, 
we can define a rotation R(6') ■ R(0) that aligns the coordinate system of the gravity-noise 
strain tensor hRf with the propagation frame of the GW. In this case, the contribution of 
the GW to the total strain tensor h assumes the simple Gartesian form 


h' — 


hx 

0 


hx 0^ 
-h+ 0 
0 0 ^ 


(7) 


In the GW propagation frame, the ; 2 -axis corresponds to the direction of propagation. The 
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two rotation matrices are given by 



^cos(6') 0 

— sin(6')^ 

R(0) = 

0 1 

0 



^sin(6*) 0 

cos(6') 

/ 


^ cos((/)) 

sin(0) 

0 ^ 

R(0) = 

— sin(0) 

cos(0) 

0 


y 0 0 ly 

and the transformation of the gravity-noise tensor reads 

= R{e) ■ R(0) ■ hRf ■ R(-0) • R{-e). (9) 

Let us now take the sum h' of a single GW and Rayleigh-wave NN ^j). 

It can be shown that the following relations hold: 

— 2cot{6)h[^ + cot^(6')h33 

-f csc^(0)27r7Gpo-^ 

( 10 ) 

hx = ^32 — COt(6*)h23 

k 

-h i csc(6')27r7Gpo^ 7 &(^) sin(ai - 0). 

i 

The sum over displacement amplitudes in the hrst equation simply denotes the total ampli¬ 
tude of vertical seismic surface displacement at the GW detector. Applying a trigonometric 
addition theorem, the sum over Rayleigh waves in the second equation can be rewritten in 
terms of horizontal seismic displacement of the Rayleigh-wave held along the two directions 
X, y. Therefore, a linear combination of tensor channels and one or two seismic channels (GM 
accelerations) can be found that perfectly cancels NN in the two target channels hn, hi 2 . 
Since the linear combination involves csc(6') and cot(0) functions that can become very large, 
it should be intuitively clear that the analytical relation cannot be used in practice when 
channels are also contaminated by additional instrumental noise. These noise contributions 
would be amplihed by the gravity-noise cancellation. A practical solution to this problem is 
investigated in Section IV. 

Next, we repeat the calculation for infrasound NN. Here, we choose to calculate the 
gravity perturbation underground. We simply want to avoid the technical problem of placing 






the GW detector inside the fluctuating density fleld (i. e. the infrasound held), which leads 
to additional terms in the NN. Avoiding these terms does not change the final results or 
applicability of the method. Below surface, the perturbation of the gravity potential by a 
plane infrasound wave reflected from the surface reads 


7Po 


( 11 ) 


where 6p{u) is the amplitude of pressure fluctuations, 7 the adiabatic coefficient of air, pq 
the mean air pressure, po fhe mean air mass density, and fch the horizontal wavenumber of 
an infrasound wave. This leads to the gravity-gradient tensor 


gis(l^= 0,uj]a,/3) = -4:7r^^6p{uj)sin^{/3) 

Wo 

^ cos^(a) cos(a) sin(a) —icos(a) 
cos(q!) sin(a) sin^(a) —isin(a) 
, —icos(a) —isin(a) —1 


\ 

/ 


( 12 ) 


The matrix is identical to the Rayleigh-wave matrix in Equation ( 6 ). As before, the angle 
a specifies the direction of propagation of the wave along the horizontal direction. The 
factor sin^(/3) is owed to the fact that the horizontal wavenumber of an infrasound wave, 
fch = A:sin(/?), depends on the angle of incidence (3 with respect to the surface normal. 

The equations for the noise cancellation are given by 


— 2 cot(6')h(3 -|- cot^(6')h33 

-f csc^(6')^— ^(5pi(a;)sin^(A), 

7Po ^ 

( 13 ) 


hx = h\2 — cot(6')h23 

-f icsc(6')^^^ sin^(/3j) sin(Q!i - 0). 

7Po ^ 

Here, we can see that the case of infrasound cancellation is more challenging. The sums 
over infrasound waves do not correspond to easily observable quantities. For example, a 
microphone collocated with the GW detector would observe independent of the 

angles a,, 0*. Directional information could come from a gravimeter sensing associated 
fluctuations of the gravity fleld. Still, the sums cannot be rewritten in terms of gravimeter 
channels (GM acceleration) due to additional factors sin(0j). Another problem of gravimeter 
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channels is that they would be dominated by seismic noise at frequencies above 10 mHz, 
which makes these channels useless for the cancellation of infrasound NN. 

In this section, we have presented analytical expressions describing a new approach to 
cancel NN in full-tensor GW detectors. Cancellation of Rayleigh NN as shown in Equation 
(10) can be achieved with tensor channels and an additional 3-axis seismometer. In the 
tensor detector described in Section II, the CM acceleration channels of the detector provide 
a three-axis linear and three-axis angular seismometer with SNR in excess of 10® at 0.1 - 
0.3 Hz. However, we have seen in Equation (13) that cancellation of infrasound NN is more 
challenging. A term remains that cannot be observed by a single microphone. Nonethe¬ 
less, it is shown in Section IV that tensor channels greatly simplify the noise cancellation. 
Furthermore, while the analytical expressions cannot be used when including instrumental 
noise, it will be shown that efficient subtraction using only a small number of environmental 
sensors is still possible, for Rayleigh and infrasound NN. 


IV. NEWTONIAN NOISE CANCELLATION IN TENSOR GW DETECTORS 

It was estimated that low-frequency GW detectors need to achieve strain sensitivities 
of about 10“^‘^Hz“^/^ above 0.1 Hz in order to have good chances to observe GWs [7]. In 
Figure 2, a sensitivity model is shown together with estimates of the seismic (Rayleigh) 
and infrasound NN. It can be seen that seismic NN needs to be suppressed by about 3 
orders of magnitude, and infrasound NN by about 5 orders of magnitude. This is a truly 
daunting challenge and is rightfully considered a potential show-stopper for ground-based, 
low-frequency detectors. In order to achieve this suppression, it was proposed that large 
arrays extending over square-kilometers made of several tens to hundreds of sensors are 
to be deployed around the GW detector. The environmental sensors need to monitor their 
signals with sufficient sensitivity to avoid signihcant sensor-noise contributions in the cleaned 
strain channels. 

Assuming array conhgurations designed with maximized efficiency (no sensor can be 
removed without increasing noise residuals to an unacceptable level), the sensor signal-to- 
noise ratio (SNR) needs to be at least as high as the inverse of the suppression goal. For 
low-frequency detectors, this means that seismometers need to sense seismic displacement 
with SNR > 1000, and microphones need to sense pressure fluctuations with SNR > 10® at 
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FIG. 2. Sensitivity target for low-frequency GW detectors as first derived for the MANGO detector 
concepts [7]. The two NN estimates are based on measured seismic [23] and infrasound spectra 

[24]. 

0.1 Hz. Requirements can be far more demanding for microphones than suggested by this 
rule of thumb. The additional challenge with infrasound NN cancellation is that the density 
perturbations are described by a 3D infrasound field, but the array can only be deployed 
at the surface. This greatly limits the ability to extract the required information from the 
infrasound measurements, and affects the optimal array configuration. Irrespective of the 
intrinsic sensitivity of microphones to pressure fluctuations, wind noise poses a challenge for 
high-sensitivity infrasound monitoring [25]. Consequently, a solution of the infrasound NN 
problem requires new methods and technology. 

It is assumed that all noise is stationary and Gaussian, which means that optimal noise 
cancellation is achieved with Wiener filters [26] . In frequency domain, the Wiener filter is a 
vector mapping reference channels R to an estimate n of the NN according to 

h{oj) = w{oj) ■ R{u). (14) 

This form makes use of the fact that noise at different frequencies is uncorrelated. If the 
Wiener hlter is applied in time domain, then the last equation needs to be substituted by 
a convolution between the filter and the reference channels [27]. The estimated NN fi is 
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subsequently subtracted from the target channel. 

A Wiener hlter is calculated from the matrix Crr of correlations between reference chan¬ 
nels, which contains the sensor noise contributions on the diagonal, and correlations Crt 
between reference channels and target channel. In general, correlations between channels 
have to be estimated from measurements, but here we assume that the density helds are 
isotropic (in a 2D sense for the Rayleigh held, and for the half-sphere of incident infrasound 
waves), which allows us to calculate the correlations precisely. Examples of calculated cor¬ 
relations between seismometers and gravity data can be found in [15, 18]. In the notation 
of the previous section, the target channels of the noise cancellation are h'i 2 - The refer¬ 
ence channels consist of all environmental sensors and the strain channels h[^, hgg. The 
components of the Wiener hlter are given by 

w{oj) = ■ (CRR(a;))“^ (15) 


In order to evaluate the performance of a Wiener hlter, we plot the residual spectrum of 
the target channel after NN subtraction relative to the initial spectrum Ctt(i^) of the target 
channel. The residual is given by [17, 18] 

r{uj) = 1 - • (<^rr(^))~^ • Ckt{u) 

Ctt(^) 

The relative subtraction residuals r{u) have frequency dependence since a distributed array 
of reference channels has frequency-dependent correlations between channels, and also since 
sensor SNRs vary with frequency, as shown in Figure 3. 

In the following subsections, we investigate Wiener hltering of infrasound and Rayleigh 
NN in detail. Results are presented only for the target channel. Noise residuals are 
similar for tint we point out that since the Rayleigh NN cancellation ideally requires 
horizontal seismometer channels, there may be additional noise from surface shear waves 
(Love waves) that contribute to horizontal surface motion without producing gravity noise. 


A. Rayleigh Newtonian noise 

In the following, we demonstrate the effect of strain reference channels on residuals after 
subtraction of Rayleigh NN. Figure 4 shows the residual noise in using 7 seismometers 
on a 5 km ring around the detector, and the vertical CM channel of the detector. The 
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FIG. 3. Estimated SNRs of reference channels. The self-noise of the seismometer lies a factor 10 
below commercially available broadband instruments. The self-noise of the microphone is already 
achieved in current instruments. Wind noise is not included (see [28] for a recent review on wind 
noise reduction methods). SOGROs GM channels function as seismometers with 1000x higher 
SNRs. 

speed of Rayleigh waves is equal to 3.5km/s assumed here to be independent of frequency. 
For simplicity, we also assume that the seismometers measure seismic displacement with 
a frequency-independent SNR = 1000 (the CM channel with 1000x higher SNR), and the 
strain channels measure Rayleigh NN with frequency-independent SNR = 1000. Frequency- 
dependent SNRs require yet to be developed numerical tools to optimize sensor arrays used 
for NN cancellation. 

It is worth discussing in detail how this result compares to the analytical expression in 
Equation (10). First, if only the seismometer at the detector were used, then residuals 
near 6 * = 0, tt would grow to values close to 1 independent of frequency. This case is 
represented in Figure 4 by the residuals at frequencies > 0.3 Hz, where the seismometers 
on the ring have vanishing impact on residual noise. Residual noise close to 1 is already 
better than predicted by Equation (10), since the cot(6*), csc(6*) factors mean that noise in 
reference channels is amplihed to inhnity at these angles. The Wiener hlter avoids excess 
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FIG. 4. Relative residuals of Rayleigh NN subtraction in channel h'^ using 7 seismometers on a 
5 km ring, and the vertical CM channel. Seismometers and strain channels have SNR = 1000. The 
CM channel is 1000 x more sensitive than the seismometers. Residuals are independent of angle (/>. 


noise as can be understood from Figure 5. It shows the non-zero Wiener hlter coefficients 
with reference channels consisting of the strain channels, the CM vertical channel, and 7 
seismometers on a 5 km ring. Since the hlter magnitude varies over orders of magnitude, 
the log-modulus transform, f[x) = sgn(a;) logio(l + |a:|), was applied [29], where sgn is the 
signum function. The dashed curves show noise residuals for inhnite sensitivity reference 
channels. As expected, hlter coefficients tend to inhnity near 6* = 0, vr. The curves diher 
from the analytical expression in Equation (10) due to additional seismometers on the ring. 
If reference channels have hnite SNRs, then the Wiener hlter has reduced coefficients towards 
the boundaries. In this way, the Wiener hlter avoids injecting excess sensor noise into the 
target channel at the price of vanishing NN cancellation. The hlter coefficients also explain 
why the noise residuals in Figure 4 are very small near 9 = ■n 12. For this value of 9, the only 
reference used in the NN cancellation is the vertical CM channel, which has very high SNR. 
In comparison, the impact of the CM channel on residuals near 6* = 0, vr is less pronounced 
since the residuals are dominated by noise from the strain and seismometer channels. 

Figure 4 also shows that, if seismometers are added to form a ring around the detector, 
then subtraction performance is increased substantially for all values of 9, especially within 
a certain frequency band. This frequency band can be chosen by adjusting the radius of 
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FIG. 5. Filter coefficients of a noiseless Wiener filter for Rayleigh NN subtraction at 0.1 Hz 
(dashed), and using reference channels with SNR = 1000 except for the CM vertical channel, 
which has SNR = 10® (solid). Coefficient w{h' 2 ^) = 0 in both cases. Gray curves in the background 
are for the 7 seismometers on a 5 km ring. 


the ring. The optimal radius for a certain target band depends on the Rayleigh-wave speed 
and on the SNR of the reference channels. The higher the SNR or the lower the speed, the 
smaller the optimal ring radius. 

A direct comparison between the case with and without tensor channels is shown in 
Figure 6. As before, the conhguration of the seismic array consists of 7 seismometers on 
a ring with 5 km radius, the vertical CM and strain channels. It can be seen that in the 
case of Rayleigh NN subtraction, strain reference channels help, but do not lead to a large 
decrease of residuals compared to the conventional scheme based on seismometers alone. 
However, one should keep in mind that the residuals based on only local reference channels 
is competitive with the residuals of the full Wiener hlter for 0.5 < 9 < 2.6. Using exclusively 
local channels not only simplihes the experimental setup, but also potentially results in an 
increased robustness of the subtraction performance with respect to wave scattering and 
contributions from local seismic sources. These claims need to be tested in more detailed 
simulations. 


15 














FIG. 6. Relative residuals of Rayleigh NN subtraction at 0.1 Hz in channel h'n using vertical 
CM and strain channels, and 7 seismometers on a 5 km ring. The solid curve shows the residuals 
including all reference channels, the dashed curve using only the 8 seismic channels, and the dotted- 
dashed curve using only local channels, which means the vertical CM and strain channels. Strain 
and seismometer channels have SNR = 1000, the CM channel SNR = 10®. 


B. Infrasound Newtonian noise 

Next, we present analogous results for infrasound NN subtraction. Noise residuals using 7 
microphones on a 1 km ring around the detector, another 7 microphones on a 600 m ring, one 
microphone located at the detector, and strain channels are shown in Figure 7. The strain 
channels measure infrasound NN with frequency-independent SNR = 10®, and we assume 
that microphones measure pressure fluctuations with frequency-independent SNR = 10^. 
With these parameter settings, noise residuals he above the required 10“® level. Increasing 
the microphone SNR and correspondingly decreasing the radii of the rings would further 
lower residuals, but developing such microphones will not be straight-forward. The combined 
effect of the two microphone rings is good broadband performance of the noise cancellation. 
Adding another smaller microphone ring would lower residuals at higher frequencies. While 
increased residuals were to be expected from Equation (13) near 6 = 0, tt, the microphones 
ensure that residuals only weakly depend on 9. Subtraction residuals at 0.1 Hz would be 
greater by almost an order of magnitude if only a single microphone ring were deployed. 
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FIG. 7. Relative residuals of infrasound NN subtraction in channel h'^i using 15 microphones. 
Microphones have SNR = 10^, and the strain channels SNR = 10^. Residuals are independent of 
angle 4>. 

Subtraction performance is generally worse compared to the Rayleigh NN case since infra¬ 
sound NN is caused by a 3D waveheld that is monitored by a microphone array constrained 
to Earth surface. Being composed of evanescent waves, the Rayleigh seismic held also dis¬ 
places the ground below surface, but since Rayleigh waves are genuinely surface waves, a 
seismic array deployed at the surface can extract all information about associated density 
huctuations. 

The Wiener hlter coefficients in Figure 8 conhrm that the infrasound NN subtraction is 
more challenging. First, the plotted hlter coefficients are far oh the analytical expression 
in Equation (13). This was to be expected since Equation (13) states that no simple com¬ 
bination of local reference channels can provide an accurate estimate of NN. In contrast 
to the Rayleigh NN case, where good subtraction can also be achieved with local channels 
only, at least for a range of values of 6, deployment of a microphone array is essential to 
achieve good subtraction performance. The dashed curves show the hlter coefficients with a 
single microphone located at the detector, while the solid curves are calculated for a Wiener 
hlter that includes 7 additional microphones on a 600 m ring. Microphones measure pressure 
huctuations with SNR = 10^, and strain channels measure infrasound NN with SNR = 10^. 
Comparison between the two sets of curves in Figure 8 shows that the microphones on the 
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ring contribute substantially to the noise cancellation, since the coefficients are greatly mod- 
ihed when including the extra microphones. Among all local channels, only hgg contributes 
to a NN estimate near 0 = 0, tt in both cases. Coefficients of the Wiener hlter for infrasound 
NN subtraction generally depend weakly on the SNR of reference channels. The coefficients 
in Figure 8 would look the same for any SNR values greater than 10. 

In order to achieve the low residuals in Figure 7, it was necessary to deploy two rings 
of microphones, a total of 15 microphones being used for the noise cancellation. At this 
point, one may wonder if the strain reference channels still contribute signihcantly to the 
subtraction performance. Figure 9 shows that, in contrast to Rayleigh NN, reference strain 
channels play a very important role in infrasound NN subtraction. The dotted-dashed 
curve shows the noise residuals if only local channels are used (i. e. strain channels and one 
microphone). Subtraction performance is very poor. The dashed curve shows subtraction 
residuals if only microphones are used. In this case, subtraction can be fairly good, but only 
for certain values of 6. If all channels are included, then excellent subtraction performance 
is achieved for all values of 6 (solid curve). This means that microphone array and strain 
channels both play an important role. 


V. DISCUSSION 

Combining Figures 6 and 9, we come to a clear conclusion towards an efficient scheme 
of NN cancellation in full-tensor GW detectors. Since Wiener hlter coefficients of the strain 
channels for Rayleigh and infrasound NN are different, it is impossible to cancel both si¬ 
multaneously. However, we have seen that Rayleigh NN can be cancelled efficiently without 
including strain channels, whereas subtraction of infrasound NN prohts greatly from strain 
channels. Consequently, we propose the following strategy. As a hrst step, Rayleigh NN 
needs to be cancelled in all 5 independent strain channels. Good broadband subtraction 
performance can be achieved with spiral seismometer arrays consisting of a few tens of seis¬ 
mometers with SNR = 1000 at the microseismic peak [7]. The required sensitivity of the 
seismometers is also near the sensitivity of available commercial instruments [30]. When the 
strain channels are cleaned from Rayleigh NN, the next step is to subtract infrasound NN 
using microphones and the strain channels. We have seen that noise residuals are greatly 
reduced compared to subtraction with microphones only. Information about the 3D infra- 
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FIG. 8. Wiener-filter coefficients using microphones with SNR = 10^ and strain channels with 
SNR = 10^ for infrasound NN subtraction. Coefficient w{h 2 ^) = 0 in all cases. Dashed: Single 
microphone located at the detector. Solid: Seven microphones on a 600 m ring around the detector 
and a single microphone located at the detector. Subtraction is calculated at 0.1 Hz. Only the 
coefficient for the microphone at the detector is plotted. 

sound field, which cannot be retrieved with a 2D microphone array is partially provided 
by the strain channels. Subtraction residuals still do not reach the sensor-noise limit, but 
compared to the scheme with microphones only, residuals are lowered by orders of magni¬ 
tude depending on the direction of propagation of the GW. Therefore, we conclude that the 
problem of NN mitigation in low-frequency GW detectors is greatly facilitated in full-tensor 
detectors. 

Gertain aspects of low-frequency NN cancellation demand or deserve a deeper analy¬ 
sis. First of all, SNRs of the sensors are frequency-dependent although we used frequency- 
independent SNRs to simplify our calculation. In a real experiment, the NN would be 
removed using Wiener hlters with sensor arrays optimized for the actual measured SNRs. 
Numerical tools for array optimization need to be developed. Second, it is conceivable that 
mitigation schemes relying on local channels rather than on data from large distributed ar¬ 
rays are more robust against wave scattering and disturbance from local sources, which both 
influence the spatial correlation of seismic or infrasound helds. Glearly, it is always easier to 
optimize a noise cancellation based on local channels, but there may be additional advan- 
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FIG. 9. Relative residuals of infrasound NN subtraction at 0.1 Hz in channel using strain chan¬ 
nels and 15 microphones. The solid curve shows the residuals including all strain and microphone 
reference channels, the dashed curve using only microphones, and the dotted-dashed curve using 
only local channels, which means the strain channels and one microphone located at the detector. 
Strain channels have SNR = 10^ and microphone channels have SNR = 10“^. 


tages of using strain channels with respect to robustness of the performance. Robustness 
can play a very important role since gravity perturbations need to be understood at a level 
1/1000 for Rayleigh NN or 1/10^ for infrasound NN. Especially for the infrasound micro¬ 
phones, it is unclear whether such level of accuracy in the monitoring itself can be achieved. 
Here, we do not refer to the intrinsic readout noise of the sensors, but environmental noise, 
for example, from wind [25]. 

Also, it is not yet understood if a change of spatial correlation due to scattering or local 
sources poses a limitation to the subtraction performance, which can only be overcome by 
deploying additional environmental sensors, or if it leads to modification of NN and density 
fields in such a way that a simple rearrangement of senors can compensate the loss in 
subtraction performance without increasing the number of sensors. These considerations 
play an important role for high subtraction goals and at low frequencies since it is highly 
unlikely that a location can be identihed where heterogeneities of the ground and surface 
profile are negligible over the relevant volumes [31]. Nonetheless, all these challenges exist 
for any type of NN mitigation at low frequencies, and one may expect that they impact 
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coherent noise mitigation more strongly in schemes based on large sensor arrays than on 
schemes partially relying on strain channels and a smaller number of environmental sensors. 

Finally, given the fact that infrasound NN needs to be understood at a level 1/10^ for a 
complete cancellation of the noise, approximations in our gravity models such as negligible 
size of the GW detector or placing the detector at the surface despite the fact that it may 
be a few hundred meters underground potentially lead to signihcant modelling errors. More 
accurate models need to be investigated to hnd out if the mitigation scheme proposed in 
this article needs to be modihed for the most ambitious noise-suppression goals. 
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